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Abstract 

We present a search for standard model Higgs boson production in association with a W boson 
in proton-antiproton collisions (pp — > W^H — > Ivbb) at a center of mass energy of 1.96 TeV. 
The search employs data collected with the CDF II detector which correspond to an integrated 
luminosity of approximately 1 fb _1 . We select events consistent with a signature of a single 
lepton (e ± /^ ± ), missing transverse energy, and two jets. Jets corresponding to bottom quarks 
are identified with a secondary vertex tagging method and a neural network filter technique. The 
observed number of events and the dijet mass distributions are consistent with the standard model 
background expectations, and we set 95% confidence level upper limits on the production cross 
section times branching ratio ranging from 3.9 to 1.3 pb for Higgs boson masses from 110 to 
150GeV/c 2 , respectively. 

PACS numbers: 13.85.Rm, 14.80.Bn 
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I. INTRODUCTION 



Standard electroweak theory predicts a single fundamental scalar particle, the Higgs 
boson, which arises as a result of spontaneous electroweak symmetry breaking [l| ; however, 
the Higgs boson has not been direct observed experimentally. The current constraint on the 
Higgs boson mass, tuh > 114.4 GeV/c 2 at 95% confidence level (C.L.), comes from direct 
searches at LEP2 experiments J2J. Global fits to electroweak measurements exclude masses 
above 144GeV/c 2 at 95% CL 

At the Tevatron pp collider at Fermilab, the next-to-leading-order (NLO) Higgs boson 
production cross section by gluon fusion is about ten times larger than for WH associated 
production, and the cross section for WH is about twice that of ZH [4|. The Higgs boson 
decay branching ratio is dominated by H — > bb for itlh < 135 GeV/c 2 and by H — > W + W~ 
for rriH > 135GeV/c 2 {{J. Background QCD bb production processes in the same invariant 
mass range have cross sections at least four orders of magnitude greater than that of Higgs 
boson production and this renders searches in the gg — > H — > bb channel extremely 
difficult. However, requiring the leptonic decay of the associated weak boson reduces the 
huge QCD background rate. As a result, WH — > tvbb is considered to be one of the most 
sensitive processes for low mass Higgs boson searches 1 . 

Searches for WH — > Ivbb at J~s = 1.96 TeV have been most recently reported by CDF 

_ h _ n 

(using data corresponding to an integrated luminosity of 319 pb 1 )[7(| and DO (440 pb )[8j. 
The CDF analysis used a secondary vertex 6-tagging algorithm (secvtx) to distinguish b- 
quark jets from light flavor or gluon jets |9j . Upper limits on the Higgs boson production rate, 
defined as the cross section times branching ratio (a ■ B), were derived for mass hypotheses 
ranging from 110 to 150GeV/c 2 . The rate was constrained to be less than 10 pb at 95% 
CL. for vtih = HO and less than 2.8 pb for 150 GeV/c 2 . In that analysis, about 50% of the 
jets tagged by the secvtx tagging algorithm were actually falsely b-tagged jets originating 
from light flavor, gluon, or charm quarks. This effect is due to the finite resolution of track 
measurements and the long lifetime of D mesons. Even the small fraction of mistagged 
events in the dominant Wqq process is significant compared to true Wbb production. To 
reduce this contamination and enhance the 6-jet purity of our sample, we introduce a b- 



1 In this paper, lepton (£) denotes electron (e ) or muon (/J,^ 1 ), and neutrino (v) denotes electron neutrino 
(e„) or muon neutrino 
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tagging neural network filter which uses as inputs jet characteristics as well as secondary 
vertex information. 

In this paper, we present a search for WH — > Eubb production at CDF using about 1 fb _1 
of data. Section [III describes the CDF II detector. The event selection criteria are explained 
in Sec. IIHI In Sec. IIV} the 6-tagging algorithm with secvtx and neural network (NN) are 
discussed in detail. Contributions from the standard model (SM) background are calculated 
in Sec. |V]for various sources. In Sec. [VT| signal acceptance and systematic uncertainties are 
estimated. The search optimization and statistical interpretation of the results are presented 
in Sees. IVIII and IVIIIl respectively. Finally, our conclusions are presented in Sec. HX1 



II. CDF II DETECTOR 



The CDF II detector geometry is described using a cylindrical coordinate system [10 ]. 
The z-axis follows the proton direction, and the polar angle 9 is usually expressed through 
the pseudorapidity rj = — ln(tan(#/2)). The detector is approximately symmetric in rj and 
in the azimuthal angle <fi. 

Charged particles are tracked by a system of silicon microstrip detectors and a large open 
cell drift chamber in the region \r]\ < 2.0 and \r)\ < 1.0, respectively. The tracking detectors 
are immersed in a 1.4 T solenoidal magnetic field aligned coaxially with the incoming beams, 
allowing measurement of charged particle momentum transverse to the beamline. 

The resolution on the transverse momentum = psin6 is measured to be Spt/pt ~ 
0.1% • pr(GeV) for the combined tracking system. The resolution on the track impact 
parameter (do), or distance from the beamline axis to the track at the track's closest approach 
in the transverse plane, is cr(do) ~ 40 fim, about 30 /im of which is due to the transverse size 
of the Tevatron interaction region. 

Outside of the tracking systems and the solenoid, segmented calorimeters with pro jective 
tower geometry are U8 ed to r ec„ n8trUCt e,ec tr0magn e tlC showers and nadronic jet s flflQ 



over the pseudo-rapidity range |r^| < 3.6. A transverse energy E? = Esin8 is measured in 
each calorimeter tower where the polar angle (9) is calculated using the measured z position 
of the event vertex and the tower location. 

Small contiguous groups of calorimeter towers with signals are identified and summed 
together into an energy cluster. Electron candidates are identified in the central electromag- 
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netic calorimeter (CEM) as isolated, mostly electromagnetic clusters which match a track in 
the pseudorapidity range < 1.1. The electron transverse energy is reconstructed from the 



electromagnetic cluster with a resolution a(E T )/E T = 13.5%/y E T /(GeV) © 2% |rjj. Jets 
are identified as a group of electromagnetic (EM) and hadronic (HAD) calorimeter clusters 
which fall within a cone of radius AR = y/ 'A<p' 2 + Ar/ 2 < 0.4 units around a high- Et seed 
cluster [l^|. Jet energies are corrected for calorimeter non-linearity, losses in the gaps be- 
tween towers, multiple primary interactions, out-of-cone losses, and inflow from underlying 
event jlol ]. 

For this analysis, muons are detected in three separate subdetectors. After at least five 
interaction lengths in the calorimeter, the muons first encounter four layers of planar drift 
chambers (CMU), capable of detecting muons with pt > 1.4GeV/c 16j. Four additional 
layers of planar drift chambers (CMP) behind another 60 cm of steel detect muons with 
Pt > 2.8 GeV/c [171 ] . These two systems cover the same central pseudorapidity region with 

< 0.6. Muons which exit the calorimeters at 0.6 < < 1.0 are tracked by the CMX 
detector, consisting of four layers of drift chambers. Muon candidates are then identified as 
isolated tracks which extrapolate to line segments or "stubs" in one of the muon subdetectors. 
A track which is linked to both CMU and CMP stubs is called a CMUP muon. 

The CDF trigger system is a three-level filter, with tracking information available even 
at the first level [l8|. Events used in this analysis have all passed the high-energy electron 
or muon trigger selection. The first stage of the central electron trigger requires a track with 
Pt > 8 GeV/c pointing to a tower with Et > 8 GeV and £'had/-E'em < 0.125. The first 
stage of the muon trigger requires a track with pt > 4 GeV/c (CMUP) or 8 GeV/c (CMX) 
pointing to a muon stub. A complete lepton reconstruction is performed online in the final 
trigger stage, where we require Et > 18 GeV/c 2 for electrons and pt > 18 GeV/c for muons. 



III. EVENT SELECTION 

The observable final state from the WH — > ivbb signal consists of two jets plus a lepton 
and missing transverse energy. The leptonic W decay requirement in WH events yields the 
high-pr lepton and large missing transverse energy due to the neutrino. 

The results presented here use data collected between February 2002 and February 2006. 
The data collected using the CEM and CMUP triggers correspond to 955 ± 57 pb _1 , while 
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the data from the CMX trigger corresponds to 941 ± 56 pb _1 . 

The missing transverse energy (Et) is a reconstructed quantity that is defined as the 
opposite of the vector sum of all calorimeter tower energy depositions projected on the 
transverse plane. It is often used as a measure of the sum of the transverse momenta of the 
particles that escape detection, most notably neutrinos. To be more readily interpretable as 
such, the raw Et vector is adjusted for corrected jet energies, for the transverse momentum 
of the muons, and for the energy deposition of any minimum ionizing high-p-r muons. 

Events are considered as WH candidates only if they have exactly one high-p T isolated 
lepton [3], with Et > 20 GeV for electrons or pt > 20GeV/c for muons. The isolation 
cone of AR = 0.4 surrounding the lepton must have less than 10% of the lepton energy. A 
primary event vertex position is calculated by fitting a subset of particle tracks which are 
consistent with having come from the beamline. The distance between this primary event 
vertex and the lepton track z Q must be less than 5 cm to ensure the lepton and the jets come 
from the same hard interaction. Some leptonic Z decays would mimic the single-lepton 
signature if a lepton is unidentified. Events are therefore rejected if a second track with 
Pt > lOGeV/c forms an invariant mass with the lepton which falls in the Z-boson mass 
window (76 < mex < 106GeV/c 2 ). The selected events are required to have Et greater 
than 20 GeV. 

The WH signal includes two jets originating from H — > bb decays; these jets are expected 
to have large transverse energy. The jets are required to be in the pseudorapidity range 
covered by the silicon detector so that secondary vertices from b decays can be reconstructed. 
Specifically, we require the jets satisfy Et > 15 GeV and \rj\ < 2.0. The search for WH — ► 
ivbb is performed in the sample of events with W+ exactly 2 jets; however, samples of events 
with W+l,3,>4 jets are used to cross-check the background modeling. 

To increase the signal purity of the W+2-jet events, at least one jet must be 6-tagged 
by the secvtx algorithm. If only one of the jets is 6-tagged, the jet must also pass the 
NN 6-tagging filter. If there are two or more secvtx 6-tagged jets, the NN is not applied. 
With a SECVTX mistag rate of 1%, it is rare that two or more jets in the same events are 
mistagged by secvtx. 
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IV. SECONDARY VERTEX b- TAGGING 



Multijet final states have dominant contributions from QCD light flavor jet production, 
but the standard model Higgs boson decays predominantly to bottom quark pairs. Correctly 
identifying the b quark jets helps to remove most of the QCD background. An algorithm has 
been developed and used to tag displaced secondary vertices from b quark decays; however, 
the sample tagged by the secvtx algorithm still has significant contamination from falsely- 



tagged light-flavor or gluon jets and the misidentification of c quarks as 6-jets |20j. This 
search introduces a multivariate NN technique intended to improve the secvtx tagging 
purity. 

The 6-quark has a relatively long lifetime, and B hadrons formed during the hadroniza- 
tion of the initial b quark can travel a significant distance on the order of millimeters before 
decaying into a collection of lighter hadrons. The decay vertex can be reconstructed by iden- 
tifying tracks which form a secondary vertex significantly displaced from the pp interaction 
point (primary vertex). 

The SECVTX 6-tagging algorithm is applied to each jet in the event, using only the tracks 
which are within r]-(j) distance of AR = 0.4 of the jet direction. Displaced tracks in jets 
are used for the secvtx reconstruction and are distinguished by a large impact parameter 
significance (l^o/^dol) where do and a do are the impact parameter and the total uncertainty 
from tracking and beam position measurements. Secondary vertices are reconstructed with 
a two-pass approach which tests for high-quality vertices in the first pass and allows lower- 
quality vertices in the second pass. In pass 1, at least three tracks are required to pass 
loose selection criteria (j>t > 0.5GeV/c, \d Q /cr do \ > 2.0), and a secondary vertex is fit 
from the selected tracks. One of the tracks used in the reconstruction is required to have 
Pt > 1-0 GeV/c. If pass 1 fails, then a vertex is sought in pass 2 from at least two tracks 
satisfying tight selection criteria {px > 1.0 GeV/c, \d /a do \ > 3.5 and one of the pass 2 tracks 
must have pr > 1.5GeV/c). If either pass is successful, the transverse distance (L xy ) from 
the primary vertex of the event is calculated along with the associated uncertainty. This 
uncertainty includes the uncertainty on the primary vertex position. Finally jets are 
tagged positively or negatively depending on the L xy significance (L xy /<TL xy ): 

Lxy/cr Lxy > 7.5 (positive tag) (1) 
L xy /a Lxy < -7.5 (negative tag) (2) 
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These values have been tuned for optimum efficiency and purity in simulated 6-jet samples 
from decays of top quarks. The energy spectrum for those jets is similar to the spectrum 
for b jets from decays of Higgs bosons. 

The sign of L xy indicates the position of the secondary vertex with respect to the primary 
vertex along the direction of the jet. If the angle between the jet axis and the vector pointing 
from the primary vertex to the secondary vertex is less than it/2, L xy is positively defined; 
otherwise, it is negative. If L xy is positive, the secondary vertex points towards the direction 
of the jet, as in true B hadron decays. For negative L xy the secondary vertex points away 
from the jet; this may happen as a result of mismeasured tracks, so jets tagged with a 
negative L xy are labeled mistagged jets. In order to reject secondary vertices due to material 
interaction, the algorithm vetoes two-track vertices found between 1.2 and 1.5 cm from the 
center of the silicon detector (the inner radius of the beampipe and the outer radius of the 
innermost silicon layer being within this range). All vertices more than 2.5 cm from the 
center are rejected. 

The negative tags are useful for evaluating the rate of false positive tags, which are defined 
"mistags" in the background estimates. Mismeasurements are expected to occur randomly; 
therefore the L xy distribution of fake tags is expected to be symmetric with respect to zero. 
Simulated events are used to correct a small asymmetry due to true long-lived particles in 
light flavor jets. 

The efficiency for identifying a secondary vertex is found to be different in the simulated 
and observed datasets. We measure an efficiency scale factor, which is defined as the ratio 
of the observed to the simulated efficiencies, to be 0.91 ± 0.06 in a sample of high-E^ jets 
enriched in b jets by requiring a soft lepton (p? > 8 GeV/c 2 ) from semileptonic heavy quark 
decays {9}. 

Secondary vertex SECVTX 6-tagging exploits the long lifetime of B hadrons. D hadrons 
originating from c-quarks also have fairly long lifetime, and secondary vertices in c-jets are 
frequently tagged. Therefore jets tagged by SECVTX are contaminated not only by falsely 
tagged light flavor (uds or gluon) jets, but also by long-lived charmed hadrons in c-jets. A 
neural network has been developed to filter the 6-tagging results in order to improve the 
6-tagging purity. 



The neural network used in this article employs the jetnet 



2l| package. The tagger is 



designed with two networks in series. The b — I network is trained to separate 6-jets from 
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light-quark jets (/-jets), and the b — c network is trained to separate 6-jets from c-jets. Jets 
which pass a cut on both of the NN outputs are accepted by the tagger. These neural 
networks are trained and applied only to events which are already tagged by the secvtx 
algorithm. The current NN 6-tagging is tuned to increase the purity of the SECVTX 6-tagged 
jets, not to increase the tagging efficiency 

The neural networks take as input the 16 variables listed in Table [B These variables 
are chosen primarily because the 6-quark jets have higher track multiplicity, larger invariant 
mass, longer lifetime and a harder fragmentation function than c- and /-quarks jets. The 
track parameters and L xy significance are good discriminators for 6-jets. The vertex p^ TX 
and invariant mass Mvtx are useful variables for identifying /-jets; however c-jets have pt 
spectra similar to fe-jets. Pseudo-cr (L xy x M VT x/pf rx ) 1 the vertex fit x 2 ; an d the track- 
based probability of a jet to come from the primary vertex are the best discriminators. The 
outputs of the two neural networks are shown in Fig. [IJ 

The NN 6-tagger is validated by comparing the performance on data and Monte Carlo 
events. The NN output from h — I network on a sample of secvtx tagged heavy-flavor 
jets from events with an electron candidate with Et > 8 GeV electron data and from the 
corresponding Monte Carlo sample are shown in Fig. [21 as are the outputs of the b — I 
network on tagged light-flavor jets from data and Monte Carlo 2 . Figure [2] shows the good 
agreement in NN 6-tagger performance between data and Monte Carlo. 

We tune the cut value for 90% b efficiency (after the secvtx efficiency), corresponding 
to a value of NN&_2 = 0.182 and NN&_ C = 0.242. The data-to- Monte-Carlo scale factor, 
measured from the electron sample, is 0.97±0.02. Note that this is an additional scale factor 
with respect to the secvtx efficiency scale factor because all of the jets under consideration 
have already been tagged by secvtx. At these cut values, the NN filter rejects 65% of 
light-flavor jets and about 50% of the c jets while keeping 90% of 6-jets after being tagged 
by SECVTX. 

2 A small but purified 6-jet sample is obtained by requiring a soft lepton in the jet. 
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SECVTX variable 



SECVTX-independent variable 



Number of tracks in fitted vertex Number of good tracks 

Vertex fit \ 2 Jet Probability [22] 

Transverse decay length (L xy ) Reconstructed mass of pass 1 tracks 

L xy significance (L xy /aL xy ) Reconstructed mass of pass 2 tracks 

Vertex Mass (M vtx = IPvtxl) 2 — (XPvtx) 2 ) Number of pass 1 tracks 

Pseudo-cr (L xy x M vtx /p^ x ) Number of pass 2 tracks 

/ (Sgood tracks Pt) Epassl track Pt/Pt 



Vertex pass number (pass 1 or 2) Ep aS s2 track Pr/p Jet 



TABLE I: Input variables used in the NN 6-tagging filter. The variables in the first column are 
properties of the identified secondary vertex, while variables in the second column are jet properties 
independent of any identified vertex. 



— bottom jets 
charm jets 
light jets 




0.4 0.6 0.8 1 
b-l Trained Network Output 



0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 




— bottom jets 
charm jets 
light jets 
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FIG. 1: Neural network outputs obtained from trainings of b vs. I jets (left) and b vs. c jets (right). 
Output distributions for b, c and I jets are shown in solid, dashed, and dotted lines, respectively. 



V. BACKGROUND 

The final state signature from WH — > iubb production can also be reached by other pro- 
duction processes. The main background processes are 1^+jets production, tt production, 
and non-W^ QCD multijet production. Several electroweak production processes also con- 
tribute with smaller background rates. In the following subsections the contribution from 
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b-l Trained Network Output b-l Trained Network Output 

FIG. 2: Comparisons of NN 6-tag output in data (solid line), and Monte Carlo (dashed line) for 
SECVTX-tagged heavy-flavor-enriched jets (left) and tagged light-flavor jets (right). 

each background source is calculated in detail. 
A. Non-iy QCD Multijet 

Events from QCD multijet production sometimes mimic the W^-boson signature with fake 
leptons or fake $t- Non-H^ leptons are reconstructed when a jet passes the lepton selection 
criteria or a heavy-flavor jet produces leptons via semileptonic decay. Non-jy $t can be 
observed via mismeasurements of energy or semileptonic decays of heavy-flavor quarks. It 
is difficult to model and produce the former class of events in detector simulation since 
the reasons for mismeasurement are not known quantitatively. Instead, we estimate the 
contribution of non-W events directly from the data sample before 6-tagging is applied. 

Generally, the bulk of non-iy events are characterized by a non-isolated lepton and small 
$t- Lepton isolation / is defined as the ratio of calorimeter energy inside a cone of AR = 0.4 
about the lepton to the lepton energy itself. The quantity I is small if the lepton is well- 
isolated from the rest of the event, as typified by a true leptonic W decay. This feature is 
used to extrapolate the expected non-VF contribution into our signal region, namely, small 
/ and large The following 4 sideband sectors are used for the extrapolation: / > 0.2 
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and $ T < 15 GeV (region A), I < 0.1 and # T < 15 GeV (region B), I > 0.2 and $ T > 20 
GeV (region C), and / < 0.1 and $t> 20 GeV (region D). Here, region D corresponds to 
the signal region. In extracting the non-H^ background contribution from data, we make 
the following two assumptions: lepton isolation and $t are uncorrelated in non-W events, 
and the 6-tagging rate is not dependent on i n non-VT events. The level at which these 
assumptions are justified determines the assigned uncertainty. 

With the first assumption, the number of non-W 7 events (Nff n ~ w ) and their relative 
fraction in the signal region before requiring 6-tagging (/ nO n-w0 obey the following relations: 

AT non-W _ N B * N C , . 

Nd ~ N A ' {6) 
_ N^- W N B xNc ( .. 

Uon-W ~ ^ NaxNd > ^ 

where Ni [i = A, B, C, D) are the number of pretag events in each sideband region. The 
number of pretag events has been corrected for known sources of prompt leptons. By in- 
voking the second assumption, the secvtx 6-tagging efficiency obtained in region B can be 
applied to the signal region D. Here we define an event tagging efficiency per taggable jet 
(one with at least two good secvtx tracks) as follows: 

^-(tagged event) 

TB = ^(taggable jet) ' ( 5 ) 

where N^ gsed evcnt - ) and j\r^ ta ss able J et ) are ^he number of tagged events and taggable jets in 
region B, respectively. Then the number of non-VT background in region D after secvtx 
b-taggmg(N^ non ~ w ) is obtained by using the "Tag Rate" relation: 

N +non-W = xrBX ^(taggable jet.) _ (g) 

It is also possible to estimate non-W contribution solely from the SECVTX-tagged sample 
as: 

AI /+non-W _ ^3 X 

~ NX ' 10 
where N£(X = A, B,C, D) in the "Tagged Method" are the number of events with positive 
tags. These methods are data-based techniques, so the estimates could also contain other 
background processes. The contributions from ti and H^+jets events to each sideband region 
are subtracted according to the calculated cross sections for those processes, including the 
appropriate tagging efficiencies. 
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To validate the four-sector method and estimate their systematic uncertainties, we vary 
the boundaries of the four regions and divide the I and Et sidebands into two E (0.1 < 
I < 0.2 and E T > 20 GeV) and F (I < 0.1 and 15<E T < 20 GeV) sidebands. The observed 
deviations imply a 25% systematic uncertainty in the non-P^ background yield, assigned 
conservatively for both the pretag and tagged estimates. 

The independent estimates from the tag rate method (Eq. [6]) and the tagged method 
(Eq. [7]) are combined using a weighted average. The result from the tagged method gives a 
slightly higher estimate than the tag rate method, but the two results are consistent within 
the 25% uncertainty. 

A non-W rejection factor associated with the NN 6-tagging filter is measured from data 
in region C. Region C has event kinematics similar to non-W events in the signal region D 
because lepton isolation is the only difference between the two regions. The non-VF estimate 
calculated before applying NN 6-tagging is scaled by this NN rejection factor; this assumes 
the NN filter is uncorrelated with the isolation. 

The non-W estimate for events with at least two secvtx tags is obtained by measuring 
the ratio of the number of events with at least one 6-tag to the number with at least two 
6-tags in sideband regions and applying the ratio to the estimate of tagged non-H^ events in 
the signal region D. 

B. Mistagged Jets 

The rate at which secvtx falsely tags light-flavor jets is derived from generic jet samples 
in varying bins of r], <p, jet Et, track multiplicity, and total event Et scalar sum. Tag rate 
probabilities are summed for all of the taggable jets in the event, jets with at least two tracks 
well measured in the silicon detector. Since the double-mistag rate is small, this sum is a 
good approximation of the single-tag event rate. Negative mistags - tags with unphysical 
negative decay length due to finite tracking resolution - are assumed to be a good estimate 
of falsely tagged jets, independent to first order of heavy flavor content in the generic jet 
sample. The systematic uncertainty on the rate is largely due to self-consistency in the 
parameterization as applied to the generic jet sample. The positive mistag rate is enhanced 
relative to the negative tag rate by light-flavor secondary vertices and material interactions 
in the silicon detectors. As a result, the positive mistag rate is corrected by multiplying 
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the negative mistag rate by a factor of 1.37 ± 0.15. This factor is measured in a control 
sample by fitting the asymmetry in the vertex mass distribntion of positive tags over negative 
tags [23]. An additional correction factor of 1.05 ± 0.03 is applied for data collected after 
December 2004, when the Tevatron beam position changed slightly. The mistag rate per jet 
is applied to events in the VT+jets sample. The total estimate is corrected for the non-iy 
QCD fraction and also the top quark contributions to the pretag sample. To estimate the 
mistag contribution in NN-tagged events, we apply the light flavor rejection power of the 
NN filter 0.35 ± 0.05 as measured using light-flavor jets from various data and simulated 
samples. 



C. W+Heavy Flavor 

The Wbb, Wcc, and Wc states are major background sources of secondary vertex tags. 
Large theoretical uncertainties exist for the overall normalization in part because current 
Monte Carlo programs generate W -(-heavy-flavor events only to leading order. Consequently, 
rates for these processes are normalized to data. The contribution from true heavy-flavor 
production in V^+jet events is determined from measurements of the heavy-flavor event 
fraction in VT+jet events and the 6-tagging efficiency for those events, as explained below. 

The fraction of VT+jets events produced with heavy-flavor jets has been studied exten- 



sively using an alpgen + HERWIG combination of Monte Carlo programs |24j . I25j . Calcula- 
tions of the heavy-flavor fraction in alpgen have been calibrated using a jet data sample, 
and measurements indicate a scaling factor of 1.5 ±0.4 is necessary to make the heavy-flavor 
production in Monte Carlo match the production in multijet data 9J. The final results of 
heavy-flavor fractions are obtained as shown in Table HT1 In the table, IB and 1C refer to the 
case in which only one of the heavy-flavor jets is detected; this happens when one jet goes 
out of the detector coverage or when two parton jets merge into the same reconstructed jet. 
Similarly, 2B and 2C refer to the case in which both of the heavy-flavor jets are observed. 

For the tagged VT+heavy flavor background estimate, the heavy-flavor fractions and 
tagging rates given in Tables [Til and II I II are multiplied by the number of pretag H^+jets 
candidate events in data, after correction for the contribution of non-iy and ti events to the 
pretag sample. 

The previous CDF analysis using 319 pb _1 of data provided some evidence that the 
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Jet Multiplicity 1 jet 



2 jets 



3 jets > 4 jets 



Wbb (IB) (%) 1.0 ± 0.3 1.4 ± 0.4 2.0 ± 0.5 2.2 ± 0.6 



Wbb (2B) (%) 



1.4 ± 0.4 2.0 ± 0.5 2.6 ± 0.7 



Wcc (1C) (%) 1.6 ± 0.4 2.4 ± 0.6 3.4 ± 0.9 3.6 ± 1.0 



Wcc (2C) (%) 



1.8 ± 0.5 2.7 ± 0.7 3.7 ± 1.0 



Wc (%) 



4.3 ± 0.9 6.0 ± 1.3 6.3 ± 1.3 6.1 ± 1.3 



TABLE II: The heavy-flavor fractions, given in percent, for the W + jets sample. The results from 
ALPGEN Monte Carlo have been scaled by the data-derived calibration factor of 1.5 ± 0.4. (Wc 
fractions have not been rescaled.) 

disagreement between the predicted and observed numbers of W+l jet and W+2 jet events is 
due to the heavy-flavor fraction [7|. In this analysis, an updated correction factor of 1.2±0.2, 
obtained by fitting tagged W+l jet events only, is applied to the heavy-flavor fraction. The 
W+ heavy flavor background contribution is obtained by the following relation: 



where fuF is the heavy-flavor fraction, e tag is the tagging efficiency, iVrop is the expected 
number of ti and single top events, and A^ewk is the expected number of WW, WZ, and Z 
boson events. 

D. Top and Electroweak Backgrounds 

Production of both single top quark and top-quark pairs contribute to the tagged lep- 
ton+jets sample. Several electroweak boson production processes also contribute. WW 
pairs can decay to a lepton, neutrino as missing energy, and two jets, one of which may be 
charm. WZ events can decay to the signal Wbb or Wcc final state. Finally, Z — > t + t~ 
events can have one leptonic r decay and one hadronic decay. The leptonic r decay gives 
rise to a lepton + missing transverse energy, while the hadronic decay yields a narrow jet of 
hadrons with a non-zero lifetime. 

The normalization of the diboson and single top backgrounds are based on the theoretical 
cross sections listed in Table [TV] the luminosity, and the acceptance and 6-tagging efficiency 



Nw+HF — /HF " e tag - [-Wpretag ' (1 — fi 



non— 
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Jet Multiplicity 1 jet 



2 jets 



3 jets > 4 jets 





> 1 SECVTX 6-tag (%) 








Whh HRl 

VV UU 1 ID 1 


9 x o | 
oo.z in z.^i 




9 


4D 9 4- 


9 Q 


Wbb (2B) 




K1 5 X 5(! KA 1 4- 

01.0 in o.u o^±.i in 


3.8 


55.1 ± 


3.9 


Wcc (1C) 




8(14-1 1 074- 


1.4 


11.6 ± 


1.6 


Wcc (2C) 




lid 4-9(1 17(1 + 


2.4 


17.8 ± 


2.5 


Wc 


Q Q _|_ 1 O 
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1.1 


3.4 ± 
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QIlr | ATAT A ton- (<%1 
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Whh HRl 

V V UU \ 1. J-J J 


29.9 ± 2.1 


31.8 ± 2.3 34.1 ± 


9 4 


35 Q 4- 


9 6 


1^66 (2B) 




47.2 ± 3.4 51.5 ± 


3.7 


51.3 ± 


3.6 


Wcc (1C) 


3.8 ± 0.5 


5.5 ± 0.8 6.1 ± 


0.9 


6.4 ± 


0.9 


Wcc (2C) 




9.9 ± 1.4 8.6 ± 


1.2 


9.5 ± 


1.4 


Wc 


5.0 ± 0.7 


4.6 ± 0.7 3.1 ± 


0.4 


3.4 ± 


0.5 



> 2 secvtx 6-tag (%) 

9.7 ± 0.7 13.6 ± 1.0 11.5 ± 0.8 
1.2 ± 0.2 0.8 ± 0.1 0.9 ± 0.1 



Wbb (2B) 
Wcc (2C) 



TABLE III: The 6-tagging efficiencies in percent for various 6-tagging strategies on individual 
W+heavy-flavor processes. Categories IB, 2B refer to number of taggable 6-jets in the events, 
with similar categories for charm jets. Those numbers include the effect of the data-to-Monte 
Carlo scale factors algorithm and the neural network filter. 



derived from Monte Carlo events 
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26 



27 



281 ] . The acceptance is corrected for lepton 



identification, trigger efficiencies, and the z vertex cut. The tagging efficiency is always 
corrected by the 6-tagging scale factor. 



E. Summary of Background Estimate 

We have described the contributions of individual background sources to the final back- 
ground estimate. The background estimates for the condition of exactly one b-tagged jet 
after applying the NN filter and at least two secvtx 6-tagged jets are summarized in Ta- 
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Theoretical Cross Sections 

WW 12.40 ± 0.80 pb 

WZ 3.96 ± 0.06 pb 

ZZ 1.58 ± 0.02 pb 

Single top s-channel 0.88 db 0.05 pb 

Single top t-channel 1.98 ± 0.08 pb 
Z r+T- 320 ± 9 pb 

tt 6.7 ±g;g pb 

TABLE IV: Theoretical cross sections and uncertainties for the electroweak and single top back- 
grounds, along with the theoretical cross section for tt at nit = 175GeV/c 2 . The cross section 
of Z° — ► r + r is obtained in the dilepton mass range m TT > 30GeV/c 2 together with a fc-factor 
(NLO/LO) of 1.4. 



bles [V] and [VI] The estimates are plotted in Figs. [3] and |U for the case of exactly one 6-tag 
before and after applying the NN 6-tag filter. The observed number of events in the data 
and the SM background expectations are consistent both before and after NN 6-tagging is 
applied. The same is true for the number of events with at least two 6-tagged jets. (See 
Table EH and Fig.H) 



VI. HIGGS BOSON SIGNAL ACCEPTANCE 

The kinematics of the SM WH — > ivbb process are well denned, and events can be 
simulated accurately by Monte Carlo programs. The pythia program was used to generate 



the signal samples 29(. Only Higgs boson masses between 110 and 150 GeV/c 2 are considered 
because this is the mass region for which the decay H —>■ bb dominates. The number of 
expected WH — > ivbb events iV is given by: 

N = e- J Cdt- a(pp -> WH) ■ B(H -> 66), (9) 

where e, / Cdt, a(pp — > WH), and B(H — > bb) are the event detection efficiency, integrated 
luminosity, production cross section, and branching ratio, respectively. The production cross 
section and branching ratio are calculated to NLO precision [5| . The acceptance e is broken 
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Jet Multiplicity 


1 iet 


2 jets 


3 jets 


> 4 jets 


Pretag Events 


94051 


14604 


2362 


646 


A/Tict q cr 


1 30 7 + 97 3 


kq q 4- 1 n 7 


157 + 31 


4.2 ± 0.8 


Wbb 


306 9 ± 106 9 


144.7 ± 49.4 


29 9 ± 9 7 


6.4 ± 2.5 


Wcc 


63 1 + 22 


43.0 ± 14.7 


8.7 ± 2.8 


1.9 ± 0.8 


Wc 


185.7 ± 47.2 


34.4 ± 9.0 


3.4 ± 0.9 


0.6 ± 0.2 


tt(6.7pb) 


6.9 ± 1.2 


42.0 ± 6.6 


84.9 ± 12.8 


98.6 ± 14.3 


Single Top 


16.7 ± 1.8 


23.5 ± 2.4 


4.8 ± 0.5 


0.8 ± 0.1 


Diboson/Z° -► t + t~ 


11.7 ± 2.2 


14.2 ± 2.3 


3.9 ± 0.9 


1.0 ± 0.3 


non-W QCD 


84.2 ± 14.1 


38.9 ± 6.7 


12.1 ± 2.3 


5.5 ± 1.2 


Total Background 


814.9 ± 140.7 


394.4 ± 66.6 


163.4 ± 18.7 


118.9 ± 14.9 


Observed Events 


856 


421 


177 


139 



TABLE V: Background estimate for events with exactly one secvtx 6-tag that passes the NN 
filter function of jet multiplicity. 



Jet Multiplicity 


2 jets 


3 jets 


> 4 jets 


Observed Events (pretag) 


14604 


2362 


646 


Mistag 


3.5 ± 0.5 


2.0 ± 0.3 


1.2 ± 0.2 


Wbb 


20.3 ± 7.0 


5.7 ± 1.8 


1.0 ± 0.4 


Wcc 


3.3 ± 1.1 


0.4 ± 0.1 


0.1 ± 0.04 


Wc 








it (6.7pb) 


10.4 ± 2.3 


29.5 ± 6.4 


45.5 ± 9.9 


Single Top 


4.2 ± 0.7 


1.4 ± 0.2 


0.3 ± 0.1 


Diboson/Z° -»• r+r" 


1.2 ± 0.3 


0.3 ± 0.1 


0.1 ± 0.1 


non-W QCD 


1.4 ± 0.3 


0.9 ± 0.2 


0.3 ± 0.1 


Total Background 


44.2 ± 8.5 


40.1 ± 6.8 


48.6 ± 10.0 


Observed Events 


39 


44 


65 



TABLE VI: Background estimate for events with at least two secvtx 6-tagged jets as a function 
of jet multiplicity. 



23 



— •— Observed 

i~ | W+Heavy Flavor 

Mi stag 
| | Non-W QCD 

Diboson/Z^Tx 

tt (6.7pb)+Single Top 
III Background Error 



W+1 jet W+2 jet 



W+3 jet W+>4 jet 
Jet Multiplicity 



|? 1200 

CD 
> 
LU 

B1000 
i_ 

Jd 

E 800 

^ i 



— •— Observed 

[_ | W+Heavy Flavor 

Mi stag 
| | Non-W QCD 
Diboson/Z°-yrc 
tf (6.7 pb)+Single Top 
Background Error 




W+1 jet W+2 jet 



W+3 jet W+>4 jet 
Jet Multiplicity 



FIG. 3: Number of events as a function of jet multiplicity for events with exactly one secvtx 6-tag 
before(left) and after(right) applying the NN 6-tagging requirement. 
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FIG. 4: Number of events as a function of jet multiplicity for events with at least two secvtx 
6-tagged jets. 



down into the following factors: 

£ — X] ( e 20 ' trigger ' Clepton ID " Cfctag " ^kinematics ' B(W — > £v)) , (10) 

where e Zo , e tr igger, eicpton id, e b tag, and e kincm atics are efficiencies to meet the requirements of 
primary vertex, trigger, lepton identification, 6-tagging, and kinematics. The major sources 
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FIG. 5: The summary of acceptance of the process WH — > iubb in W+2jet bin for various 6-tagging 
strategies as a function of Higgs boson mass. 

of inefficiency are the lepton identification, jet kinematics, and 5-tagging factors; each is 
a factor between 0.3 and 0.45. The factor e Zo is obtained from data, and the others are 
calculated using Monte Carlo samples. The total signal acceptances for various 6-tagging 
options including all systematic uncertainties as a function of Higgs boson mass are shown 
in Fig. [5l 

The expected number of signal events is estimated by Eq. [9] at each Higgs boson mass 
point. The expectations for various 6-tagging strategies are shown in Table IVHl The NN 
6-tagging filter keeps about 90% of signal while it removes 35% of the total background in 
W+2 jet events as shown in Fig. [3l 

The total systematic uncertainty on the acceptance stems from the jet energy scale, ini- 
tial and final state radiation, lepton identification, trigger efficiencies, and (^tagging. A 2% 
uncertainty on the lepton identification efficiency is assigned for each lepton type (CEM elec- 
tron, CMUP and CMX muon), based on studies of Z boson events. For each of the high pt 
lepton triggers, a 1% uncertainty is measured from backup trigger paths or Z boson events. 
The initial and final state radiation systematic uncertainties are estimated by changing the 
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parameters related to ISR and FSR from nominal values to half or double the nominal 
The difference from the nominal acceptance is taken as the systematic uncertainty. The 
uncertainty in the incoming parton energies relies on the eigenvalue uncertainties provided 
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Higgs Mass 


Expected Signal Events 




(GeV/c 2 ) 




1 tncr with NNtno 


> 2 tae 


110 




1 SV _|_ n 1 8 
l.o / in W.lo 


66 ± 13 


115 


Q QQ_|_n OR 1 QH X fl 1 ^ 


1.00 It U.IO 


54 + 11 


120 


3.23±0.23 1.45 ± 0.12 


1.26 ± 0.12 


44 ± 09 


130 


2.05±0.15 0.93 ± 0.08 


0.81 ± 0.08 


0.28 ± 0.06 


140 


1.03±0.07 0.46 ± 0.04 


0.40 ± 0.04 


0.15 ± 0.03 


150 


0.40±0.03 0.18 ± 0.02 


0.16 ± 0.02 


0.06 ± 0.01 



TABLE VII: Expected number of WH — ► £z/66 signal events with systematic uncertainties for 
various 6-tagging options, where "tag" and "NNtag" stand for secvtx 6-tagging and NN 6-tagging, 
respectively. 

in the PDF fits. An NLO version of the PDFs, CTEQ6M , provides a 90% confidence interval 
of each eigenvector [3l| . The nominal PDF value is reweighted to the 90% confidence level 
value, and the corresponding reweighted acceptance is computed. The differences between 
nominal and reweighted acceptances are added in quadrature, and the total is assigned as 
the systematic uncertainty 9J. 



The uncertainty due to the jet energy scale uncertainty ( JES) [15] is calculated by shifting 
jet energies in WH Monte Carlo samples by ±lcr. The deviation from the nominal accep- 
tance is taken as the systematic uncertainty. The systematic uncertainty on the secvtx 
6-tagging efficiency is based on the scale factor uncertainty discussed in Sec. IIVI When NN 
6-tagging is applied, the scale factor uncertainty is added to that of secvtx in quadra- 
ture. The total systematic uncertainties for various 6-tagging options are summarized in 
Table IVnH 



VII. OPTIMIZATION OF SEARCH STRATEGIES 

The search strategy is optimized by calculating a signal significance defined as S/ \^B, 
where S and B are the number of expected signal and background events. In this analysis, 
S and B are counted within a window which gives the best significance in dijet mass dis- 
tribution for the particular Higgs mass hypothesis being considered. The window itself is 
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source 



uncertainty (%) 



1 Tag 1 Tag & NNtag > 2 Tag 



1 ot™\'1" r~\~f\ 1 1 ) 


9 n 


/o 


-£.U /o 


9 n% 


irigger 




/o 


<^ J- /o 


^1/0 






/U 


1 8% 


4 3% 

■-± . O / o 


FSR 


2.8 


% 


3.2% 


8.6% 


PDF 


1.6 


% 


1.7% 


2.0% 


JES 


2.3 


% 


2.3% 


3.0% 


6-tagging 


3.8 


% 


5.3% 


16% 


Total 


5.8 


% 


7.2% 


19% 



TABLE VIII: Systematic uncertainties for various 6-tagging requirements. The labels "Tag" and 
"NNtag" refer to secvtx and NN 6-tagging, respectively. 

optimized by varying the window peak and width for each 6-tagging strategy. A comparison 
of the significance for various 6-tagging options, shown in Fig. [6l provides an a priori metric 
that predicts which selection gives the best result. 

Requiring the NN filter improves the sensitivity by about 10% in the sample of events 
with exactly one b tag. The significance in double-tagged events is almost the same as 
that in events with at least one tag and no NN filter. Combining the two results therefore 
yields another sensitivity improvement. This combined use of two separate 6-tagged samples 
provides a significant improvement as shown in Fig. El The total significance increases by 
20% moving from "> 1 tag" to separate categories "1 tag w/ NNTag" and "> 2 Tag." 
Therefore, we quote final results from events having exactly one secvtx 6-tagged jet passing 
the neural network filter or at least two secvtx 6-tagged jets. 

VIII. LIMIT ON HIGGS BOSON PRODUCTION RATE 

As shown in section IVHt there is no significant excess number of events over the SM 
background expectation. Because the dijet mass resonance is a useful discriminant for the 
Higgs boson signature, we use a binned likelihood technique to fit the observed dijet mass 
distributions in Figs. [7] and El and set an upper limit on the WH production cross section 
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FIG. 6: Comparison of significance obtained from various 6-tagging strategies. "Tag" and "NN 
Tag" represent secvtx and NN 6-tagging, respectively. The filled circles correspond to the com- 
bined analysis which treats disjoint samples with exactly one NN 6-tag and two secvtx tags 
separately. 

times H —>■ bb branching ratio. 

A. Binned Likelihood Technique 

The number of events in each bin follows the Poisson distribution: 

Pi(n t ,fa) = ^-— (z = l,2,---,iV bin ) (11) 

where rii, fa, and iVbin represent the number of observed events in the i-th bin, the expectation 
in the i-th bin, and the total number of bins. The Higgs production hypothesis is constructed 
by setting fa to fa = Si+bi, where Sj and 6j are the number of signal and expected background 
events in the i-th bin. This quantity Sj can also be written as a product: 

s t = a(pp - W ± H) ■ B(H -> bb) ■ e WH ■ j Cdt ■ ff H (12) 

where / 4 is the fraction of the total signal which lies in the i-th bin. In this case, a(pp — > 
W ± H) ■ B(H — > bb) is the variable to be extracted from data. An upper limit on the Higgs 
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FIG. 7: Dijet mass distribution in W+2 jets events including exactly one secvtx 6-tagged jet that 
passes the NN rj-tagging filter. The contributions of the various background sources are shown 
in histogram, while the hatched box on the background histogram represents the background 
uncertainty. 

boson production cross section times branching ratio a(pp — ► W ± H)-B(H — > bb) is extracted 
by using a Bayesian procedure with a likelihood denned by: 

l= n^.^)= n^V (13) 

i=l i=l n i- 

The background prediction 6j includes contributions from the various background sources 
described in Sec. [V] 

ft. = N TOP fJ OP + N® CD ff CD , (14) 

where ff OP and f® CD are the fractions of the total number of top (including ti and single 
top) and QCD backgrounds (including W+jets, non-W, and diboson) in mass bin i. There 
are systematic uncertainties in the estimates of both the number of signal events and the 

expected background. Such uncertainties modify the likelihood to be 

m - 



cr 



B ) = / / / n 



nQ° d Jn top Jn wh i ml 
i=i 1 



x G(N QCD , a QCD )G(N TOP , a TOP )G(N WH , a WH )dN QCD dN TOP dN WH (15) 

where the G(N, cr) factors are truncated Gaussian densities constraints using the estimated 
numbers of events and the associated uncertainties. We assume a uniform prior for a ■ B 

29 



o 

> 12 

CD 

O 

8 10 



CD 
> 
LU 



-»- Observed 
| | W+Heavy Flavor 
| | Mi stag 
| | Non-W QCD 
B Diboson/Z°-> tc 
■ B (6.7pb)+Single Top 

Background Uncertainty 
— WHx10 (m H =115GeV/c 2 ) 




50 100150 200 250 300 350 400 450 500 
Dijet Mass (GeV/c 2 ) 



FIG. 8: Dijet mass distribution in W+2 jets events including at least two secvtx 6-tagged jets. 



and integrate the likelihood over all parameters except a • B. A 95% credibility level upper 
limit on a ■ B is obtained by calculating the 95* /l percentile of the resulting distributions. 

To measure the expected sensitivity for this analysis, background-only pseudo- 
experiments are used to calculate an expected limit in the absence of Higgs boson production. 
Pseudo-data are generated by fluctuating the individual background estimates within total 
uncertainties. The expected limit is derived from the pseudo-data using Eq. [151 

The likelihoods from events with exactly one secvtx 6-tagged jet passing the NN b- 
tagging filter and events with at least two secvtx 6-tagged jets criteria are multiplied 
together. The systematic uncertainties associated with the pretag acceptance, luminosity 
uncertainty, and uncertainty of the 6-tagging efficiency scale factor are considered to be 
100% correlated between the two selection channels. Background uncertainties, specifically 
on the heavy-flavor fractions and 6-tagging scale factor, are also completely correlated. The 
"=1 tag w/ NNtag" selection combined with ">2 Tag" gives the best expected limit, as 
expected from the sensitivity study (see Fig. [6]). 

The observed limits as a function of the Higgs boson mass are shown in Fig. [9] and Ta- 
ble IPX} together with the expected limits determined from pseudo-experiments. An ensemble 
of limits from pseudo-experiments and the observed limit for each Higgs boson mass point 
are shown in Fig. [TD1 The limit in the low mass region is at most two standard deviations 
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FIG. 9: 95% confidence level upper limit on a(pp — > WH)-B(H — > bb) with an integrated luminosity 
of 1 fb _1 obtained from the combined likelihood between events with exactly one secvtx 6-tag 
passing the NN b-tagging and events with at least two secvtx 6-tagged jets. The previous CDF 
data 7] and recent DO data [^] are shown for comparison. 



higher than the expected limit, but this is consistent with a statistical fluctuation in the 
dijet mass distributions (see Fig. [7j) around m# = 115GeV/c 2 . Such a fluctuation is much 
larger than the expectation for SM Higgs boson production in this channel. 

The search sensitivity is improved significantly with respect to previous searches, about 
30% beyond the expectations from simple luminosity scaling. The two main effects are the 
separation of the 6-tagged data sample into single- and double-tagged events, and the NN 
filter applied to the single-tag sample. 



IX. CONCLUSIONS 

We have presented a search for the standard model Higgs boson in the tvbb final state 
expected from WH production. The event selection includes an additional neural network 
fe-tag filter to reduce the background contributions from light flavor and charm quark jets. 
This improvement, along with a total dataset corresponding to 1 fb _1 , allows us to improve 
the upper limit on Higgs boson production. We set a 95% confidence level upper limit on the 
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FIG. 10: Results of 95% confidence level limits obtained from the combined likelihood in pseudo- 
experiments. The arrows indicate the observed limits. 
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fTiffEfs Mass 
GeV/c 2 


Upper Limit (pb) 
Observed Expected 


SM 
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o o_i_n q 


0.16 


115 


3.4 


2.2±0.8 
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120 


2.5 


2.0±0.7 


0.10 


130 


1.6 


1.8±0.7 
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1.4 


1.7±0.6 
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150 


1.3 


1.5±0.6 


0.011 



TABLE IX: Observed and expected upper limits on a{pp -► WH) ■ B(H -► 66) at 95 % C.L., 
compared to the SM production rate calculated at NNLO. 



production cross section times branching ratio varying from 3.9 to 1.3 pb for Higgs boson 
masses 110 to 150GeV/c 2 . 
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